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FIELD OF THE INVENTION 
This invention relates to color printing, and more particularly, to a process for 
designing optimized halftone screens formed of combinations of non-orthogonal 
supercells. 

BACKGROUND OF THE INVENTION 

Halftoning or screening is the reproduction of grayscale images using dots of a 
single shade but varying size to simulate the different shades of gray. The pattern of dots 
in a halftone screen may be coarse or fine. The frequency of a screen is usually known as 
its screen ruling and is generally defined in "lines per inch (lpi)". A screen ruling of "60 
lpi" means that sixty dots lie along a linear inch, measured along the axis which brings 
the dots closest together, its screen angle. When a digital marking engine is used to 
produce a screen pattern, the available grid of machine pixels is divided into tiles called 
"halftone cells", each of which cell contains a number of machine pixels. A pattern of 
hafltone dots is created by getting the marking engine to mark some of the machine 
pixels within each halftone cell, leaving others unmarked. 

Each halftone cell is typically a two dimensional array of predetermined threshold 
levels. The halftone cells are spatially replicated and tiled to form a halftone screen. 
Halftone screening compares the required continuous color tone level of each pixel for 
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each color separation with one of the predetermined threshold levels. If the required 
color tone level is darker than the threshold halftone level, a color spot is printed at the 
specified pixel. Otherwise the color spot is not printed. Typically, the shape and tiling 
geometry of the halftone cell is a square, rectangle, parallelogram, line, or the like. 
5 Because of the constraints of the machine pixel grid, not all arbitrary rotations of 

halftone screens can be implemented digitally. There are many methods of generating 
rotated halftone screens. In the rational tangent method, each halftone cell is exactly the 
same size and shape as its neighbors and contains the same number of machine pixels. 
However, there are relatively few combinations of screen angles and frequency that work 
H 10 ^d the quantization errors of angle and frequency can be quite large. 
% It is important to have a large number of possible screen angles to select from in 

f|l full color printing using overlaid halftone screens for each of the cyan, magenta, yellow 

m and black separations. In addition to quantization errors, in two, three and four color 

printing visible moire may result from the interference between overlapping dot patterns. 
'« 15 Also, different types of images frequently require different types of halftone screens. 
q Since a particular halftone screen may not work for all image types in a particular output 

*M job, much effort has gone into the design of optimized halftone screens. For example, 

p various digital halftone screens having different shapes and angles are described in "An 

Optimum Algorithm for Halftone Generation for Displays and Hard Copies," by T. M. 
20 Holladay, Proc. Soc. for Information Display, 21, p. 185 (1980). 

Supercells have been used to increase the angular accuracy of the halftone screen 
and to reduce the number of tone levels that the screen can produce. A supercell is 
generally an array of halftone cells, however not all of the individual halftone cells in the 
supercell need be the same size or shape. The supercell fits over the machine pixel grid 
25 using the rational tangent method. 

Halftone screens formed of supercells offer designers a greater number of screen 
angles from which to choose for finding a particular optimized solution. However, 
searching for supercells can be time consuming and current techniques are limited to 
searching orthogonal-shape supercells, therefore, do not provide a complete set of 
30 solutions. It is desirable to have more design techniques which provide an even greater 
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number of screen angles. It is desirable to have a technique for finding supercells which 
can be in general parallelogram shape. 

SUMMARY OF THE INVENTION 

Halftone screens formed of non-orthogonal supercells, where both the supercell 
and the subcells are defined by arbitrary spatial vectors, provide the halftone designer 
with a new source of optimized halftone screens. A method of constructing a halftone 
screen includes selecting a frequency and screen angle of interest. A subcell having 
spatial vectors which satisfy the selected frequency and screen angle of interest is 
identified. A supercell comprising an array of the subcells is formed. An integer 
relationship potentially having numerous solutions exists between the supercell and the 
subcells. The integer relationship is solved for values of the integers and then tested 
against the values for the subcell spatial vectors. Although the solution may in some 
cases be the null set, in many cases there will be numerous solutions. Each resulting 
solution, if any, is then tested according to any additional constraints or tolerances 
specified for the particular halftone screen. If any of the resulting supercell solutions 
satisfies the tests, that supercell may be used to create a halftone screen. 

Supercells may be generated for any desired frequency and screen angle 
constraints and for any other constraints that can be translated into a set of spatial vectors. 
If the desired constraints pertain to moire-free color halftoning, for example, the number 
of solutions using non-orthogonal supercells is frequently an order of magnitude greater 
than the number of solutions using orthogonal supercells. 

DESCRIPTION OF THE DRAWINGS 
Figure 1 is two-dimensional spatial vector representation of a single-cell halftone 

cell; 

Figure 2 is a frequency vector representation of the halftone cell of Fig. 1; 
Figure 3 is a two-dimensional spatial representation of a supercell hafltone screen; 
Figure 4 is a two-dimensional spatial representation of two tiled supercell halftone 
screens; 



Figures 5 is a flow chart outlining one exemplary embodiment of a method for 
constructing a supercell halftone screen; 

Figure 6 is flow chart outlining a method of halftoning and rendering a color 
image; and 

5 Figure 7 is a block diagram of an image processor operative to perform the 

method of Figure 6. 

DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS 
As shown in Fig. 1, a single-cell halftone screen can be an arbitrarily shaped 
10 parallelogram, which can be represented in the spatial domain by two vectors, Vi = (xi, 
yi) and V2 = (X2, y2>. In other words, the two spatial vectors Vi and V2 are specified by 
spatial coordinate values xi, yi, X2 and y2. If the given single-cell halftone screen 
represented by the graph shown in Fig. 1 is used by a halftone screening system, such as 
that shown in Fig. 7, the output will appear as a two-dimensional repeated or tiled pattern. 
15 It should be appreciated that, by using Fourier analysis, a two-dimensional 

repeated pattern can be represented in the frequency domain by two frequency vectors, Fi 
Ly = (fxi, fyi) and F 2 = (fx2, frc). Fig. 2 illustrates, for example, the two frequency vectors, 

Fi and F2, corresponding to the Fourier transforms of the two spatial vectors, \% and V2, 
respectively. It should also be appreciated that the two frequency vectors Fi and F2 are 
20 perpendicular to the two spatial vectors Vi and V2, respectively, and that their vector 
lengths are given by: 

N = ^d< and (la) 
1 (ib) 



m 



1 Z| l^sin^l 

where 8 represents the angle between the vectors Vi and V2, which is also equal to the 
25 angle between the frequency vectors Fi and F 2 . In general, the moduli IFil and IF 2 I are 
real numbers. 

Also, the area A of the parallelogram encompassed by V2 and Vi can be 
represented as: 

A = [VjV 2 sin0| (2a) 
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Alternately, the area A of the parallelogram can be written as a function of the 
spatial coordinates Xi, yi, x 2 and y 2 ; i.e., as: 

A = \x 1 y 2 -x 2 y 1 \ (2b) 

Correspondingly, using Eqs. (2a) and (2b), Eqs. (la) and (lb) can be recast as: 
Ivi' 



1^1 =— , and (3a) 



A 



N = — • (3b) 

A 



The vectors Fi and F 2 can be decomposed into their scalar components as: 

ft fxi= ^ L - (4a) 

I f ^ = A' (4b) 

10 f x =^2., and (4c) 

\1 2 A 



Substitute xi, y lf x 2 and y 2 in Eq. (2b) by Eqs. (4a) - (4d), the area A of the 
parallelogram can be also written in terms of frequency components, i.e., 

A = 



/^fxlfyl fxlfyl ' ^ ^ 



15 Therefore, Eqs. (4a) - (4d) express the frequency-to-spatial-component 

relationship for a cell defined by the spatial vectors Vi and V 2 . The frequency 
components, f x i, f y i, f^, and f y2 , can be completely defined by the four coordinate values, 
xi, yi, x 2 and y 2 , and vice versa. Since Eqs. (4a) - (4d) describe a corresponding 
"mapping" of the frequency components to the spatial components, it should be 

20 appreciated that any analysis of the desired conditions in the frequency domain can be 
easily translated into a spatial domain specification. It should be appreciated that, while 
the above equations are developed in relation to a non-orthogonal single-cell halftone cell 
having a parallelogram-like shape, it is considered that orthogonal parallelograms, for 
example, squares, rectangles, etc, are special cases of general non-orthogonal 
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parallelograms and above equations are certainly suitable to describe orthogonal 
parallelogram shaped cells. 

Referring to Figs. 3 and 4, a supercell is defined by two arbitrary spatial vectors 
Ui = (mi, nO and u 2 = (m 2 , n 2 ) that form a parallelogram, where mi, n x , m 2 and n 2 are 
integers. The supercell is an array of subcells, which are defined by two vectors vi = (xi, 
yi) and v 2 = (x 2 , y 2 ) that form the parallelogram where xi, yi, x 2 and y 2 are real numbers. 
A method described herein finds all possible supercell parallelograms that have the 
property that when the supercell parallelogram is tiled, the subcell can also be tiled. 
Figure 4 shows two supercells defined by vectors ui and u 2 which are tiled. 

Both the supercell and the subcell are generally represented by parallelograms. It 
should be noted, however, that both the supercell and the subcell can be any arbitrarily 
shaped polygon that is capable of being tiled in a periodic array. The arbitrarily shaped 
subcell can be fitted within the parallelogram. 

A method of constructing non-orthogonal supercells includes selecting a subcell 
based on a desired frequency and screen angle and then finding all the supercells that can 
be formed using the desired subcell. Finding all supercell solutions includes searching a 
generalized relationship between the supercell and subcell for integer solutions. Once a 
series of integer solutions are determined, the solutions are tested against the original 
requirements and any other specified constraints. Supercells can be designed for any 
desired properties that can be transformed into spatial vector relationship, including 
moire-free conditions. The methods described herein may be used in combination with 
the methods described in co-pending, co-assigned U.S. patent application number 
09/698,104 filed October 30, 2000, U.S. patent application number 09/909,319 filed July 
13, 2001 and U.S. patent application serial number 09/905,556 filed July 13, 2001. 

The inventors found that a simple relationship exists between the subcell and the 
supercell: 

kivi + k 2 v 2 = ui 5(a) 

k 3 vi + tjv 2 = u 2 5(b) 
where ki, k 2 , k 3 and k* are integers. Unfortunately, as of the time of the invention, 
equations 5(a) and 5(b) cannot be solved analytically. It is possible that, in the future, 
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advanced computing and analytical resources may be developed that will enable an 
analytical solution to equations 5(a) and 5(b). Furthermore, equations 5(a) and 5(b) 
frequently do not have any solutions. Based on current analytical tools and computing 
power available, most solutions involve approximating Vi and v 2 . Thus for most practical 
applications, based on the current state of technology, a tolerance level is usually given as 
part of the design parameters. 

For the supercell in Figs. 3 and 4, the corresponding equations are: 

Ui = 5vi - 2v 2 and 

u 2 = 2vi + 3v 2 

Equations 5(a) and 5(b) can be converted into scalar equations: 
mi = kixi + k2x 2 6(a) 
ni = kiyi + k 2 y 2 6(b) 
m 2 = k 3 xi + kpc 2 6(c) 
n 2 = k 3 yi +k4y 2 6(d) 

Referring to Figures 5, equations 5(a) and 5(b) can be (approximately) solved by a 
direct searching. Direct searching involves specifying values for certain parameters of 
the equations, then solving for the remaining parameters. 

In step 100 desired frequencies and screen angles are specified by two frequency 
vectors, Fj = (f X i, fyi) and F 2 = (fjQ, fy2). The corresponding spatial specification by 
vectors vi and v 2 is obtained by solving Eqs. (4a) - (4d) in step 110. 

In step 120 we select a set of integer values, ki, k 2 , k 3 and Lt, which will be 
incremented to the maximum value K, respectively. Note that ki to can be positive 
and negative integers and we loop through all possible ki to combinations. For a given 
set of ki, k 2 , k 3 and k4, we calculate Eqs. (5a) and (5b) in step 130 to find the exactly 
desired supercell solution, ui and u 2 , which in general are real-number specified. Next, 
in step 140, an approximate integer solution of the supercell is found by rounding off the 
real-number vectors, Ui and u 2 , to the closest integer vectors, u'i(mi, nO and u 2 (m 2 , n 2 ), 
where mi, m, m 2 and n 2 are integers. To check if the approximate solution, u'i and u 2 , 
satisfies the tolerance requirement, in step 150, the corresponding subcell, v\ and v 2 is 
calculated by solving Eqs. (5a) and (5b) with integer vectors u'i and u 2 . In step 160, the 



approximate solution v'i and v 2 is compared with the desired subcell Vi and v 2 . If the 
difference is within the tolerance, the supercell solution, u'i and u 2, is saved (step 170), 
otherwise, the searching loop continues with another set of integers, ki, k 2 , k 3 and (step 
180). 

The foregoing generated supercells can be used in a halftone screen in a method 
of halftoning and rendering. Referring to Fig. 6, in summary, a method 1910 of 
halftoning and rendering a color image, that solves a predetermined requirement, can 
begin with an analization step 1912. In the analization step 1912 various requirements, 
such as process effects, characteristic of a target rendering device, particular image type 
requirements or other requirements, are considered. For example, process effects such as 
directions and frequency components of effects, such as, for example, development order 
effects, dual beam effects, photoreceptor velocity non-uniformity effects, mirror wobble 
effects, and/or raster start position jitter effects, may be determined. 

In a screen set selection step 1914 a search is carried out to find a set of halftone 
screens having fundamental screen frequencies and screen angles. The search may be 
carried out using the method described with reference to Figures 5A - 5D in order to find 
a set of halftone screens formed of supercells. 

Once a screen set has been selected, colorants or color separations are associated 
with individual screens from the selected set, in screen association step 1916. When 
colorants have been associated with respective screens from the selected screen set, an 
image is received in a color image reception step 1918. The image may be received from 
any image source. For example, the image may be received directly from an image 
authoring tool, via a computer network, from digital media, such as a magnetic or optical 
disk, or from the output of a document scanner. Whatever the source, the image is either 
in a raster format or a rasterized version is created through known rasterization processes. 
In the rasterization process or subsequent to the rasterization process, the image is 
separated into color planes related to the colorants previously associated with screens 
from the selected screen set. 

In a rendering step 1920, each color separation is halftoned through the halftone 
screen of the related colorant and rendered in a rendering step 1922. The image is 
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rendered through the use of a rendering device, such as, for example, a xerographic 
printer, a lithographic printer, an inkjet printer or other printer or display device. For 
example, the method is beneficially applied to at least some image-on-image xerographic 
printers and in many lithographic environments. 
5 Referring to Fig. 7, an exemplary image processor 2010 operative to carry out a 

method 1910 of halftoning and rendering a color image in accordance with a 
predetermined set of optimization parameters, includes a screen set repository 2014, a 
halftoner 2018, an image input device 2022, a temporary image storage device 2026, a 
long term or bulk storage device 2030, general image processor components and 
■H 10 functions 2034, a rendering device 2038 and a screen set searcher 2042. 
sQ The image input device 2022 can be any image data source. For example, the 

LlJ image input device 2022 can be a scanner, electronic media reader such as, a disk drive, 

fll cdrom drive, or computer communications network connection. Whatever the source, 

m image data is stored in the temporary image storage device 2026. The temporary image 

f 15 storage device 2026 is, for example, a computer memory. 

O The screen set repository 2014 includes a non-volatile memory device. The 

J?l screen set repository 2014 stores and provides access to a selected set of halftone screens, 

U such as, a set of screens selected in the screen set selection step 1914. The selected 

screens are predetermined and loaded into the screen set repository. The screen set 
20 selection may be based on characteristics and idiosyncrasies of, for example, the print 
engine 2038 or the screen selection may be based on predetermined optimization 
requirements, such as image type. Alternatively the screen set may be selected in real 
time by the screen set searcher 2042. The screen set searcher 2042 selects a screen set 
based, for example, on rendering preferences determined, by a processor operator or 
25 other expert or agent (not shown). Again, the selection may be based on the 
characteristics of the print engine 2038 or some other defined criteria. The halftoner 
2018 reviews the image data in the temporary storage device 2026 and communicates 
with the selected screen set repository 2014 in order to compare the image data with 
threshold information stored in the selected screen set. The results of that comparison are 
30 passed, for example, to the general image processor functions 2034 and form the basis for 
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the generation of control signals for operating the print engine 2038. Alternatively, the 
results of the comparison are stored in the bulk storage device 2030 for later processing. 

The general image processor functions 2034 are known in the art to include 
editing and image enhancement functions as well as print engine communication 
functions. The general image processor 2034 may modify the information delivered from 
the halftoner 2018 or simply pass it directly to the print engine 2038. 

The print engine can be any image-rendering device. In a xerographic 
environment, the print engine is a xerographic printer, such as an image-on-image 
xerographic printer. In a lithographic environment the print engine may be a lithographic 
printer, including, for example, printing plates etched through the use of the selected 
screen set. Xerographic printers are known to include a fiiser, a developer and an 
imaging member or photoreceptor. Alternatively, the print engine 2038 may be another 
device, such as an electronic display or an inkjet printer. 

Of course, image processors that are operative to carry out the method 1910 of 
halftoning can be implemented in a number of ways. In the exemplary image processor 
2010, halftoner 2018, and the general document processor functions 2034 are 
implemented in software that is stored in a computer memory, and run on a 
microprocessor, digital signal processor or other computational device. Other 
components of the document processor are known in the art to include both hardware and 
software components. Obviously the functions of these modules can be distributed over 
other functional blocks and organized differently and still represent an embodiment of the 
invention. 

The invention has been described with reference to a particular embodiment. 
Modifications and alterations will occur to others upon reading and understanding this 
specification taken together with the drawings. The embodiments are but examples, and 
various alternatives, modifications, variations or improvements may be made by those 
skilled in the art from this teaching which are intended to be encompassed by the 
following claims. 
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